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We consider a model in which massive stars form in a self-gravitating accretion disk around an 
active galactic nucleus (AGN). These stars may evolve and collapse to form compact objects on 
a time scale shorter than the accretion time, thus producing an important family of sources for 
LISA. Assuming the compact object formation/inspiral rate is proportional to the steady-state gas 
accretion rate, we use the intrinsic hard X-ray AGN luminosity function to estimate expected event 
rates and signal strengths. We find that these sources will produce a continuous low-frequency 
(< 1 mHz) background detectable by LISA if more than ~ 1% of the accreted matter is in the form 
of compact objects. For compact objects with masses > WMq, the last stages of the inspiral events 
should be resolvable above a few mHz, with rates as high as a few hundred per year. 
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I. INTRODUCTION AND MOTIVATION 

Some time in the near future, the launch and operation of the Laser Interferometer Space Antenna (LISA) l] 
will issue in a new age of gravitational wave (GW) astronomy. One of the major challenges that will face GW 
astronomers is the successful discrimination between instrumental noise, stochastic GW backgrounds, and individual, 
resolvable GW sources. For example, the population of galactic white dwarf binaries in close orbits will provide a 
major contribution to the LISA noise curve in the 0.1 — 1 mHz band Q. At the same time, this confusion "noise" can 
also be treated as a signal, and its shape and amplitude will provide important information about the distribution 
and properties of white dwarf binaries in the galaxy (e.g. [3j). 

It has recently been proposed that, in the self-gravitating outer regions of accretion disks of active galactic nuclei 
(AGN), massive stars could form and evolve, eventually collapsing into compact objects and merging with the central 
black hole 0, |]| . This final inspiral stage would be an important source of GWs in the LISA band. Depending on the 
specific properties of the accretion disk and the embedded compact objects, this population could contribute to the 
GW background and also produce individual resolvable signals. It is therefore a matter of theoretical and practical 
interest to understand the nature of such a population. 

In this paper we attempt to derive a relationship between the observable electromagnetic (EM) emission and the 
predicted GW emission from AGN out to cosmological distances. In particular, we use the hard X-ray luminosity 
function of Ucda et al. 6] to infer the accretion history of supermassive black holes (SMBHs) out to redshifts of 
z 3. Then we assume a few simple scaling factors, such as the average (Eddington-scaled) accretion rate and the 
efficiency of converting accretion energy to X-rays, and derive the time-averaged GW spectrum that might be seen by 
LISA. Another important model parameter is the black hole spin, which will determine the radiative efficiency for gas 
accretion as well as the GW efficiency for inspirals. Based on relativistic MHD simulations of accretion disks 0,11,111, 
as well as some recent observations [lfj , there is growing evidence that most AGN black holes should be rapidly (but 
not maximally) spinning, with a/M ~ 0.9 — 0.95. 

Depending on the specific model parameters, we find this background could be an important class of LISA sources, 
similar in strength and event rates to extreme mass-ratio inspirals (EMRIs) from captured compact objects ^lj. As 
in those sources, here too it is a matter of preference as to whether the steady-state background should be thought 
of as signal or noise. But for higher masses (perhaps as large as m ~ 10 5 M Q Q, provided such objects are not 
tidally disrupted before producing significant GW power, disk-embedded compact objects should produce individual, 
resolvable inspiral events with high signal-to-noise over a wide band of frequencies. 

The outline of the paper is as follows: in Sec. [H] we give a brief overview of notation, and describe the cosmological 
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X-ray luminosity distribution from AGN. In Sec. IHII we discuss the parameterization of the GW signal, from a single 
inspiral event to an integral of all sources over redshift and AGN luminosity. We then present the main results, 
showing the predicted GW power spectra for a range of model parameters and also for parameters that may depend 
on SMBH mass and redshift. In Sec. II VI we estimate the expected event rates and duty cycle for LISA, which will 
be used to help distinguish between stochastic signals and resolvable ones. In Sec. we discuss implications for the 
LISA mission as well as future GW observatories. 



II. THE HARD X-RAY LUMINOSITY FUNCTION OF ACTIVE GALACTIC NUCLEI 



We begin with a short discussion of notation. The results derived below include a number of dimensionless parame- 
ters, most of which can take values between and 1. We divide these parameters into three general classes: efficiencies, 
fractions, and densities. Efficiencies, denoted by rj, are believed to be determined by more basic physics, and typically 
have more stringent lower- and upper-limits. Fractions, denoted by /, are more model-dependent parameters and 
less-well known than the efficiency parameters, and thus have a larger range of acceptable values. Lastly, various 
cosmological density parameters fl are given as fractions of the critical density p c = 3Hq/(8ttG). They are most likely 
constrained to very small values relative to unity. A summary of all the dimensionless parameters used in the paper 
appears in Table [IJ along with acceptable and preferred values. 

The rest mass in supermassive black holes is estimated as fisMBH — 2 x 10 -6 , determined from galaxy surveys using 
the scaling relation between velocity dispersion and SMBH mass (the well-known "M-cr" relation; [LllllL lllL The 
energy density in the diffuse X-ray background between 2 keV and 10 keV is estimated as fix - 2.6 x 10~ 9 [lfij. to 
which AGN are believed to make a considerable contribution Q. Now, assume that a fraction / acc < 1 of the SMBH 
rest mass density is due to accreted gas of cosmic density f2 acc (as opposed to mass gained through mergers), which 
releases electromagnetic radiation with an efficiency of rj enl < 1. We further assume that a fraction fx < 1 of the 
electromagnetic output is in the form of X-rays in the range 2 — 10 keV. This implies 

/acc ^SMBH — (1 — ffcm) ^acc (1) 

and 



^X — ( — ) /unobsc T?X ^acc , (2) 

\l + z/ 

where rjx = fxVcm < 1, and the average inverse redshift of the sources ((1 + z) -1 ) — 0.4 is due to the redshifting of 
radiation energy. The parameter / un obsc in the above equation is the fraction of the intrinsic X— ray luminosity that 
escapes the AGN region without obscuration or reprocessing. Since from observations Ox/Osmbh - 1.3 x 10~ 3 , this 
gives the constraint 

1 \ /unobsc /acc ??X ^ ^X ~ ' 1 3 x 10 3 (3) 



, 1 + Z I 1 - 77 cm SlsMBH 

Since, for standard accretion disk theory rj CTa < 0.3 [l6[ . this immediately implies 

/unobsc /acc fx Vcm > 2 X 10~ 3 (4) 

and thus 

2 X 10~ 3 < /unobsc, /acc, ??X, /x, Vcm < 1 ■ (5) 

The lower limit only depends on the fact that AGN contribute a considerable fraction of the observed X-ray back- 
ground. 

A growing consensus has been forming that SMBHs grow almost exclusively by accretion, determined by linking 
luminosity distributions to the SMBH mass distribution, suggesting / acc ~ 1. lia llJ, lia llH |20l |2J| . A corollary of 
this assumption is that most AGN should be rapidly spinning, with r) em ~ 0.15 — 0.3 for dimensionless spin parameters 
of a/M ~ 0.9 — 0.998 16]. As mentioned in the introduction, recent MHD simulations suggest an upper limit to 
the spin parameter, due to magnetic torques that remove angular momentum from the inner edge of the disk, thus 
preventing the accreted matter from spinning up the black hole to maximal spin 0, 0, • Thus we set the fiducial 
spin parameter at a/M = 0.95, corresponding to rj cm ~ 0.2. 

For /unobsc — 0.3 6] , Eq. © would then require an X-ray efficiency of accretion of rjx — 7 x 10~ 3 . This would require 
either that the total radiative efficiency rj em would have to be of the same order, or that most of the electromagnetic 
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TABLE I: Glossary of dimensionless parameters, with allowable and preferred values 



symbol 


min 


max 


preferred 


description 


face 





1 


1 


fraction of SMBH mass due to accreted gas 


fco 





1 


0.01 


fraction of accreted matter in form of compact objects 


.fx 





1 


0.03 


fraction of EM radiation in X-rays 


/Edd 





>1 


0.1 


typical fraction of Eddington luminosity /accretion rate 


/unobsc 





1 


0.3 


fraction of emitted X-rays not absorbed/reprocessed 


?]cm 





1 


0.2 


accretion efficiency of converting gas to EM radiation 


TfK 







0.006 


accretion efficiency of converting gas to X-rays 


Igw 





l 


0.2 


accretion efficiency of converting compact objects to GW 










radiation 


^acc 





^SMBH 


2 x 10 _b 


fraction of critical density in accreted gas 


fix 





firad 


2.6 x 10~ 9 


fraction of critical density in X-rays (2-10 keV) 







Hm 


2 x 10~ 6 


fraction of critical density in SMBHs 







^SMBH 


2 x 10 -10 


fraction of critical density in GWs 


Qm 





1 


0.3 


fraction of critical density in matter 


firad 





1 


5 x icr 6 


fraction of critical density in radiation 







1 


0.7 


fraction of critical density in vaccuum energy 



emission is emitted in other bands. While such low radiative efficiencies are certainly possible (e.g., ADAF models 
of Ref. during the periods of largest growth and thus highest AGN activity, the accretion disks should be 

radiatively efficient. Indeed, the X-ray fraction fx is governed by bolometric corrections and is of the order of 0.03 if 
AGN emission is dominated by infrared and optical frequencies 0, |2jJ, Hj| . This is consistent with the fact that 
the background energy in the X-ray band from 2 to 10 keV is dominated by AGNs, whereas the energy density in the 
infrared background which is about a factor 200 higher, is dominated by ordinary galaxies. 

Pending a fuller understanding of the star formation mechanism in the accretion disk, for now we simply assume 
that a certain fraction / co < 1 of the accreted material is in the form of compact objects which will not get tidally 
disrupted before plunging into the SMBH. Since the astrophysical parameters that actually determine this fraction 
are not well known, we set it to a conservative value of 0.01. If it were much higher, the disk would be entirely 
fragmented and thus not efficiently emit EM radiation. And as we will see below, a value much below 1% would result 
in a GW signal undetectable by LISA. We further assume that a certain fraction n gw of the rest mass m of these 
compact objects is emitted in GWs during the inspiral event. Neglecting the plunge and ringdown stages (as well as 
magnetic torques in the innermost disk), we will generally set ry gw = rj em . 

Following the above notation, we now write the X-ray luminosity Lx as a fraction fx of the bolometric luminosity, 
which in turn is a fraction /Edd of the Eddington luminosity LEdd : 

Lx = /x/Edd^Edd(M) = /x?7cmM acc c 2 , (6) 

where the gas accretion rate is M a cc- The Eddington limit is a function only of the SMBH mass: LEdd(-ftl) = 
1.3 x 10 38 (M/M Q ) erg/s. Over the range of redshifts and luminosities we are probing, typical accretion rates during 
the period of maximum black hole growth are estimated to be /Edd ~ 0.1, but could likely be even greater than unity 
Ellll El . This rate is also derived from a comparison of the luminosity distribution and the mass distribution, via 
the efficiency parameter rj crn . Note that Eq. 10, together with fx ~ 0.03 and /Edd ~ 0.1, implies that this luminosity 
function corresponds to SMBH masses 1O 5 M < M < 10 10 Mq. This is consistent with typical SMBH masses inferred 
from velocity dispersion observations. Note, however, that /Edd varies during the lifetime of an AGN and, strictly 
speaking, is a distribution for each SMBH mass and redshift. In the present work, we will understand /Edd as an 
average during the active periods of AGNs, when they actually contribute to the X-ray luminosity function used 
below. 

Emission in the infrared and hard X-rays is less obscured than in other bands and are thus more easy to observe. 
Since infrared luminosity functions are more poorly known, we will use the X-ray luminosity function to obtain a 
realistic picture of AGN distributions. We stress that all our luminosities are to be understood as intrinsic, i.e. as 
a measure of the X-rays pro duced directly by the accreting gas, and before reprocessing and/or partial absorption 
within the host galaxy [^.l20j| . This is most natural for our purposes as we are dealing with electromagnetic and GW 
emission by accretion, so we are concerned primarily with emission from the very inner-most regions of the disk. 

The luminosity function is defined as the comoving number density n of objects per logarithmic luminosity interval, 
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with units of Mpc 3 . We parametrize the intrinsic X-ray luminosity function per comoving volume after Ref. 

A 



Lx-j^-(z,Lx) = 
atix 



(l + z) Pl 



z < z c (L x ) , 
z > z c (L x ) , 



(7) 



^ [1 + Z C (L X )] P 2-P2 

for L m [ n ~ 10 415 ergss -1 < Lx < 10 46 5 ergss -1 ~ L max . For a given Lx, the cut-off redshift is 

z c {Lx) - 




(8) 



10 44 ergs s 1 , 



7i = 0.86, 



and the best fit parameters in Eqs. Q and JSJ) are given by A = 2.2 x 10 -6 Mpc -3 , L* 

72 = 2.23, px = 4.23, P2 = —1.5, z* = 1.9, L a = 10 ergss -1 , and (3 — 0.335 (note that our value for the parameter 
A differs from Ref. [f| by a factor of (In 10) due to our slightly different definitions of the luminosity function). The 
distribution function defined by Eq. (7J) is plotted in Fig. Q] for a range of redshifts. The "anti-hierarchical growth" of 
AGN is evident from the increasing average luminosity with redshift: at early times, larger black holes were growing 
fastest, and most smaller black holes (M ~ 10 6 — 10 7 M©) are thus relatively young |20ll2ll|. 




L x (erg/s) 



FIG. 1: Analytic fits to the hard X-ray luminosity function given in Ref. 
thus AGN mass) increases with redshift. 



[. Note that the average intrinsic luminosity (and 



With this luminosity function, we can calculate directly the diffuse energy density in hard X-rays due to AGN: 

dn 



POO 


dt 




\ dz 




1 + zJ 


10 


dz 





dL x L x 



dU 



(9) 



In Eq. 10, cosmology enters through \dt/dz\ = [(1 + z)H(z)] 1 and, for a flat geometry, 



h(z) = h [n M (i + zf + n A ] 



1/2 



(10) 



Throughout this paper we will assume a flat, ACDM Universe with Om = 0.3, 17a = 0.7 (all other contributions to the 
energy density are negligible), and Hq = 72 km s _1 Mpc -1 |27|. Typically, we will integrate the luminosity function 
out to redshifts z max ~ 3. 

Reference has demonstrated that, when integrated over redshift and correcting for absorption and reprocessing, 
the luminosity function (J7J| explains practically all of the diffuse extragalactic X-ray flux between ~ 1 keV and a few 
hundred keV measured by the HEAOl, ASCA, and Chandra observatories. Combining Eqs. Q and Eq. (JSJ, we obtain 
fix — 8.7 x 10 _9 / uno bsc, implying / un obsc — 0.3. This is consistent with the estimate that about 60% of AGNs of 
luminosity Lx ^ 10 44 ergss -1 are obscured [2(|- The unobscured fraction of these AGN dominates fix- 

More information may come from future X-ray missions such as Astro-E2/HXD, NeXT, Constellation- X, and 
XEUS. With their larger collecting area and greater sensitivity, these instruments will allow us to extend the X-ray 
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luminosity function to lower luminosities and higher redshifts. Future IR observatories like the James Webb Space 
Telescope will greatly improve our estimates for fx and funohsa: There may also be significant AGN emission in soft 
7— rays, as recently reported from INTEGRAL observations [23]. Future 7— ray missions such as GLAST will further 
help us understand the AGN energy budget at this end of the spectrum. 



III. THE GRAVITATIONAL WAVE SIGNAL 



A. Individual inspiral spectrum 

We consider the inspiral of a single compact object of mass m into a SMBH of mass M ^> m, motivated by scenarios 
such as those described in the introduction 0,O| . It is not known whether inspiraling black holes in accreting discs 
can acquire considerable eccentricity by interacting with the disc or with other orbiting masses, and if this eccentricity 
will have time to decrease due to GW radiation reaction as the small body spiral in. Considering the astrophysical 
uncertainties in deriving the GW spectrum, we limit the analysis to the simplest inspiraling model, i.e. the circular 
equatorial one. 

Using geometrized units such that G — c — 1, a geodesic particle on an equatorial, circular orbit around a Kerr 
black hole has an orbital frequency (as measured by an observer at infinity) of |29| 



2 7T (r A ' 2 ± aVM) 



and specific energy 

E(r) 



r(r 2 -3Mr±2aVMr) 1 /2 



(12) 



We estimate the total energy emitted in gravitational waves down to a radius r as E gw (r) = m — E(r). The energy 
emitted in GWs between frequency / and f + df for such an event is 



dE gw dE gw ( df 



df dr \ dr 



(13) 



We restrict the GW emission to the leading quadrupole formula, thus we consider only GW frequencies twice the 
orbital frequencies. In the Newtonian limit, we reproduce the well-known result (e.g. [30ll3l|) valid for circular orbits 

f^(f) = J W) 2/3 - (14) 

We will generally want to restrict Eq. I|13|) to a range of frequencies / m i n < / < / max . Here, / m ; n = is the 
smallest resolution frequency-bin determined by the mission mission lifetime T Q b s and / max is the GW frequency at 
the inner-most stable circular orbit (ISCO). The location of the ISCO is given by the radius where dE gV! /df in Eq. IT5|> 
vanishes. It strongly depends on the spin of the black hole, giving / max ~ 4 — 20 mHz for a SMBH mass 10 6 M Q [29j . 
Note that rj gw = df(dE gw /df)/m is given by Eq. (JT3J|. 

In Fig. Elwe show the GW spectrum for a single inspiral event for a range of black hole masses and spins, plotting 
the so-called characteristic GW amplitude h c (f) for single events, defined by 

(compare with Figs. 3-7 in Ref. |32|), where r(t) is the comoving coordinate, dr = (1 + z)dt. 



B. Integrated time-averaged spectrum 



To estimate the total contribution from all sources, we start by integrating the GW energy density (per logarithmic 
frequency) over redshift z, averaged over time scales large compared to all physical time scales related to the sources. 
Following Ref. j^H], we have 



dp g w(f) 
din/ 



dz 



l + z 



dE gw ( f 
dfz 



(16) 
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FIG. 2: Characteristic GW strain amplitudes for individual inspiral events, where a black hole with m — IOMq merges with 
a SMBH of mass M = 10 6 , 10 7 , 1O 8 M0 at a redshift of z = 1. For each value of M, we show the spectra for two spin values, 
a/M = (solid) and 0.95 (dashed). The dot-dashed and dotted lines are the sky-averaged LISA noise curves with and without 
the contributions from galactic binaries, respectively. 



where f z = (1 + z)f, and R(z) is the rate of inspiral events per comoving volume. Note the similarity to Eq. ©, the 
energy density in hard X-rays. 

To derive an expression for R(z), we observe that the event rate for a single AGN is simply the X-ray luminosity 
divided by the total X-ray energy emitted between inspiral events and thus for a given luminosity Lx, the event rate 
is 

Assuming equal efficiencies rj gw = rj cm , this X-ray energy can be written as 

fx 
t 

J CO 

Combining Eqs. (|17f) . (|18f) and integrating over the luminosity distribution function, we get 



Ex — E gw — . (18) 





fco 


POO 


dt 


d\nf 


fx 


Jo 


dz 



fx J d\nL x E gw 

Combining with Eq. Ijlfcifl and integrating over redshift, the total (time-averaged) gravitational wave spectrum is 

rr- 1^ ^irrT-wT^' (20) 

The GW spectrum E gw (f) from each individual AGN is a function of the SMBH mass, which in turn is determined 
by the X-ray luminosity through Eq. JJj}. Note that the integrated spectrum is independent of m, as long as m is 
small enough so that the inspiral waveforms cannot be individually resolved. One measure of this resolvability is the 
duty cycle, described in the next section. 

Following fiefs. 0, 112, we will want to compare directly the stochastic background defined in Eq. I|20|l to the 
spectral density of the detector noise S n (f), which has units of inverse frequency. In this case, \/ fS n (f ) will have 
units of dimcnsionless strain. Averaging over the entire sky, weighted by the LISA antenna pattern, gives 

4 1 dpll(f) 
tt f 3 din/ 

Throu ghou t the paper we use the so-called sky and detector averaged instrumental spectral density for LISA, given 
by [H13: 

S instI (f) = (6.12 x 1CT 51 f- 4 + 1.06 x 1(T 40 + 6.12 x 1(T 37 / 2 ) Hz" 1 , (22) 
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augmented by the white-dwarf galactic confusion noise 



, \-7/3 

5f 1 (/) = 1.4xlO- 44 ( T ^) Hz^ 1 , (23) 



and the white dwarfs extra-galactic confusion noise 



S£-**(f) = 2.8 x 10" 46 (j^j ' Hz_1 • ( 24 ) 
Thus, the total (instrumental plus confusion) noise is 

S n (f) = min {s instr (/)/exp (-KT^ SBio JN/df) , S^if) + Sf(f)} + S^ a \f) . (25) 
Here, dN/df is the number density of galactic white-dwarf binaries per unit gravitational- wave frequency given by 

dN o i /l Hz\ 11/3 , . 

w = 2 x 10- Hz- . (26) 

C. Dependence on model parameters 

In this section we show the impact of each of the model parameters on the predicted GW spectrum. In this approach, 
all parameters are taken as constant in time and independent of the other parameters and the AGN masses. In Fig. [3 
we show the time-averaged GW spectrum calculated from Eq. 12(11 , for the fiducial model parameters listed in Table 
H] (solid black curves), along with the effects of varying the individual parameters around these baseline values. 

In the top left panel of Fig. [3] we show the dependence on /Edd relative to the fiducial model. For smaller values of 
/Edd) the observed luminosity function will imply a SMBH mass distribution function shifted towards higher masses 
and thus lower frequencies (blue curve) . Conversely, higher values of /Edd require smaller black hole masses and thus 
higher frequencies (red curve). 

The top right panel of Fig. [3] shows the dependence on / co , the fraction of accreted matter in the form of compact 
objects. The linear dependence of f2 gw (/) on / co is reasonable: more compact objects give more signal, which is 
also clear from Eq. (|20J) . Similarly, the predicted GW power increases when the fraction fx decreases, as seen in the 
bottom left panel of Fig. [3] The largest GW signal in this scenario is achieved for minimal fx- This corresponds to 
only a small fraction of the total accretion activity in the Universe being observable in hard X-rays. Thus we infer 
that the bolometric luminosities are actually much higher, which in turn implies more matter being accreted and thus 
more compact objects for a fixed f co . The fraction fx also determines the AGN mass cutoff through Eq. ©, since 
(fixing all other parameters) a smaller fx implies a larger Eddington luminosity, and thus higher mass and lower 
frequency. 

Motivated by Refs. 0,0], we have used a large (yet not maximal) spin parameter of a/M — 0.95 for most of the 
calculations here. The bottom right panel of Fig. shows the fiducial model along with the corresponding signal 
for a/M = 0, 0.7, and 1.0. Despite the significant difference between the signals from individual inspiral events (see 
Fig-El, the difference in total, time-averaged power is rather smaller. This is because rapidly spinning black holes 
are more efficient (high ?7 cm and ?7gw), so for a given luminosity, there is less gas being accreted and therefore fewer 
embedded compact objects, and less GW power. Similar to the results in the bottom left panel of Fig. [21 varying 
a/M also changes the inferred SMBH mass for a given luminosity, in turn changing the cutoff frequency. 

We should mention that the GW signal appears to depend insignificantly on the maximal redshift z max for z max > 3, 
where the luminosity function becomes more and more uncertain. Another parameter that we have explored is the 
low-luminosity cutoff of the luminosity function Eq. JJJ). By extending the cutoff to lower luminosities, we include 
lower-mass AGN, and thus higher frequencies. After trying L min = 10 40 5 ergs s _1 and = 10 39 5 ergs s _1 we 

found that, as expected, the lower luminosities give more signal at higher frequencies, but not significantly within the 
LISA band. However, this may be an important factor when designing a GW observatory with more sensitivity in 
the ~ 0.1 — 10 Hz band. On the other hand, the smaller accretion disks and shorter time scales of low- mass AGN 
may make them insignificant sources of disk-embedded compact objects, and thus the somewhat artificial cutoff of 
L m in = 10 41 ' 5 ergs s _1 may actually be physically justified. 
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FIG. 3: Time- averaged GW spectra for the baseline model parameters listed in Table U] (solid black curves), and the effects 
of changing the various parameters: /Edd (top left panel), f co (top right), fx (bottom left), and a/M (bottom right). The 
dot-dashed and dotted lines are the sky-averaged LISA noise curves with and without the contributions from galactic binaries, 
respectively. 



D. Sensitivity to varying model parameters 

Up to now we have assumed that the fiducial values of the parameters in Tableware constants independent of SMBH 
mass and/or redshift. Particular parameters are, however, likely to vary throughout evolution. There is evidence, for 
example, that the Eddington ratio f^dd was higher in the past than it is now, see, e.g., Refs. [U|3i|. In order to 
assess the possible influence of such evolution, following Ref. [3i| wn0 studied an AGN sample from the Sloan Digital 
Sky Survey, we model the dependence of the average f^dd on the SMBH mass and redshift as follows (while restricting 
0.01 < / E dd < 1.0): 

where j(M) is given by Eq. (4), Table 1, and Fig. 4 in Ref. |3^|. This is, in fact, probably more realistic than 
the constant /Edd — 0.1 approximation, because the SMBH density resulting from the X-ray luminosity function 
at small redshift z ~ 0.01, fisMBH — 1-3 x 10 -6 is closer to the one inferred from the velocity dispersions in the 
AGNs: Note that ^smbh ~ / dLx(dn/dLx)M/p c , where M is related to Lx via the generalization of Eq. © to 
varying / Edd , L x = fxhdd{M, z)L Ed d(M), and M is restricted to 10 5 M Q < M < 10 10 M©. The time-averaged GW 
spectrum resulting from Eq. I|27(l is shown in the left panel of Fig.0] Compared to the fiducital model (black curve), 
when varying /Edd (red curve), the slope becomes smaller, i.e. the signal at low frequency decreases whereas at high 
frequencies it increases, consistent with the tendencies seen in the top left panel of Fig. [3] Overall, the spectrum 
seems relatively robust with respect to variation of /Edd- 

Furthermore, since the larger mass disks seem more susceptible to fragmentation into smaller mass objects, / co 
might increase with M , whereas m might decrease with M. For large M , the longer accretion times also make it more 
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likely for lower-mass stars to form and evolve to compact objects before reaching the inner edge of the disk. On the 
other hand, the feedback into the disk from star formation may increase gravitational stability and thus decrease f co . 
For purpose of illustration we consider each of the following scalings: 



M 



/ co (M) = 1CT J {joSM^ ) and (28a) 

/ co (M) = 10- 3 f^^V' 5 , (28b) 



V M 

in all cases requiring / co < 0.5 to maintain a coherent gas accretion disk in order to produce the necessary X-ray flux. 
The resulting time-averaged GW spectrum is shown in the right panel of Fig. When f co increases with M, the 
signal is considerably enhanced at low frequencies (blue curve) , while the opposite occurs when f co decreases with M 
(red curve). 




FIG. 4: Time- averaged GW spectra for the baseline model parameters listed in Table U] (solid black curves), along with models 
where /Edd varies according to Eq. 1271 (left panel, red curve) or the fraction / co of accreted matter in form of compact objects 
varies according to Eq. I|28|l (right panel, red and blue curves). 



IV. EVENT RATES, DUTY CYCLE AND CONFUSION NOISE 

The total inspiral rate as seen from Earth at frequency / can be written as (see, e.g. Ref. |3l|') 

J Q l + zdz /x1c m m Jo (1 + z)H(z) J Lm . n d\nL x 

where dV/dz = Airr 2 (z) / H (z) is the fractional volume element, the Hubble rate at redshift z is given by Eq. i|10[l . 
and r(z) is the comoving coordinate, dr — (1 + z)dt. In the second expression of Eq. 1)29(1 we have used Eq. 1(19(1 with 
£" gw = r] gw m and i m ax(/) is the maximum luminosity for which the associated SMBH mass emits to frequencies up 
to /. Note that for a fixed compact object fraction / co , the rate Eq. (|29|) is inversely proportional to the typical mass 
m of the inspiraling compact object. 

The duty cycle T>(f) at a given frequency is the average number of sources contributing at any given time. This 
can be estimated by multiplying the integrand of Eq. Ij29(l with the time scale over which a source radiates around the 
frequency / with a coherent phase development. In the case of adiabatic, circular inspirals, the system emits GWs at 
a well-defined frequency f(t) which evolves monotonically in time. In the local rest frame (i.e. ignoring cosmological 
redshifts) we can estimate the coherence time as [34| : 

, m _ / 5 f W e M Q \ 2/3 ( 10*M Q \ ( 10- 3 Kz \ 8/3 
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The duty cycle as observed at Earth will thus be estimated by multiplying the integrand of Eq. 129fl with (l+z)t co h(/z) : 

= r ^pil dz f dL ^*i_ tcoh{fz) . (31) 

JxVomm J H{z) J GflnLx 

Note that it is proportional to f C o/(fyJ1cmjn 2 ), with the extra power of m coming from the dependence of i co h in 
Eq. (J20J). Smaller m means more compact objects, and also slower inspiral rates, thus each source spends more time 
around a given /. 




FIG. 5: Duty cycle £>(/) for the nominal parameter values and a range of compact object masses m. When T>(f) < 1, the 
inspiral signals should be individually resolvable. The high frequency cut-off is due to the somewhat arbitrary cut-off of the 
X-ray luminosity function at low luminosities which depends on the experimental sensitivity. 

In Fig. |S] we show the duty cycle for the fiducial model parameters and a range of compact object masses m. The 
cutoff around 20 mHz is due to the somewhat artificial low-end cutoff in the luminosity function, corresponding to 
a minimum value for M and thus maximum attainable frequency. When the duty cycle is < 1, a detector will see a 
non-gaussian, non-continuous signal which rather has a "popcorn" character ffi. Provided the signal to noise ratio 
(SNR) is sufficiently large, the detector will then be able to resolve individual events of duration t co h(f), occurring at 
a rate T(f) < t^ h (f). 

When the duty cycle is > 1, several sources are present at any given time at a frequency / in a bandwidth A/ ~ /. 
However, the number of sources in the smallest resolution frequency-bin, which corresponds to 1/T b s ~ 10~ 8 Hz, is 
not larger than one. Thus, our background is in principle subtractable if the SNR is sufficiently high. 

Let us estimate, as done in Ref. |ll|, the unsubtractable portion which will constitute the confusion noise S coni (f) 
[or dp™ i {f)/df related to S'f^ ni {f) by the analogue of Eq. (|21|) ]. For each inspiral source observed over the LISA 
mission lifetime, we evaluate the SNR defined by 

/"/max U2( f\ 

SNR 2 = / dlnf^l, (32) 

where h c (f) is the characteristic strain defined in Eq. (|15JI and S n is given by Eq. I|25|l. This integral is performed 
over the frequency range / m i n < / < /max through which the system sweeps during the mission lifetime T Q b s - While 
in principle the inspiraling object may reach /isco & n d then plunge sometime in the middle of the LISA mission, for 
concreteness we only consider systems that are active throughout T b s - In Fig.[S]we show SNRs for a few typical cases 
as a function of the time t remaining to plunge (as of the end of the observation). We then insert a step- function 
factor 6 (SNRthr - SNR) in the integrand of Eq. (23) to obtain dp c ™ l (f)/df. This assures that only events with SNR 
too small to be subtracted out contribute to this confusion noise. 

We note in passing that, strictly speaking, unresolvable events contribute to the S n noise and thus the SNR in 
Eq. I|32(l should be obtained by an iterative procedure in which S n includes events that could not be subtracted out. 
This becomes important when the merger noise of events whose SNR is smaller than the threshold SNR is comparable 
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FIG. 6: Signal to noise ratio for the inspiral of a compact object of masses (left panel) IOMq and (right panel) lOOM© onto a 
SMBH of mass M = 10 5 M (solid line), M = 10 6 M© (dashed line), and M = 10 7 M & (dot-dashed line) at 1 Gpc distance, as 
a function of the time left to the ISCO at the end of the observation time T b s = 3yr. Also shown is the threshold of SNR=15 
(dotted line). 



or larger than the instrumental noise. In this case we will slightly underestimate the fraction of the true unsubtractable 
background. 

With the SNR calculation in hand, we can also estimate the inspiral rate above a given threshold. In the left panel 
of Fig. [7] we show both the total inspiral rate in the Universe (solid curves), as well as the rate of detectable signals 
for SNR > 15 (dashed curves), for to — 1 — 10 3 M Q . The total rate, which here can be thought of as the rate of 
compact objects reaching the ISCO, is given at each redshift as 

rw - (33 , 

1 + z az 

We have also defined [see Eq. 1)29(1 ] an inspiral rate as a function of frequency which is a measure of the number 

of chirping systems that passes through a frequency / per unit time. To distinguish these rates, we will refer to T(z) 
as the "plunge rate" and as the "chirp rate." 

While the plunge rate is proportional to m _1 , the SNR increases with to, so the smaller mass inspirals are only 
resolvable out to smaller distances. Also, at a given redshift, the integrated SNR generally decreases with increasing 
M as the lower /isco resides in the region of higher instrumental noise. By integrating the area under the dashed 
curves in Fig. [7] we find the resolvable event rates to be ~ 150yr _1 for to = 1 — 1OM and ~ 50yr _1 for to = 100Af Q . 
Inspirals with to — 10 3 M Q , while observable out to high redshifts, would be relatively rare, with rates ~ 15yr _1 . 
From Eq. (|29|) . we see that these plunge rates are proportional to f co , so could potentially be used to constrain that 
poorly-known parameter. 

In the right panel of Fig. we plot the chirp rate T(f) for the same range of compact object masses. At low 
frequencies, the chirp rate funtion is nearly flat, due to a steady-state "flux conservation" as each inspiraling object 
enters a frequency bin, another will leave it. Then, at higher frequencies, systems begin to drop out altogether as 
they reach the plunge frequency for a given mass M, up to the final cut-off frequency due to the lower limit end of 
the SMBH mass function (compare with the duty cycle T>(f) in Fig. [SJ. 

Note that these event rates are, within orders of magnitude, consistent with independent estimates of EMRIs based 
on loss cone calculations. Such estimates give inspiral rates of "captured" compact objects per galaxy ranging from 
~ 10~ 7 j3|| to ~ 10 _5 yr _1 37]. With a few billion galaxies out to z ~ 1, this corresponds to rates between ~ 100 
and ~ 10 4 yr _1 . 

In Fig. [H] we show the unresolvable portion of the GW power spectrum for the fiducial model and a range of to, 
subtracting out all contributions from systems with SNR^15. Clearly, for larger values of to, the threshold SNR is 
reached at an earlier time and thus lower frequency in the inspiral, reducing the confusion noise at higher frequencies. 
We notice that even the subtracted events can increase the total noise l|25[l by an amount proportional to S n [38|. 
because information used to subtract the events cannot be used to detect other events. However, in our case, we have 
found that this increase in negligible. 

From the dependence T>(f) cx fco/m 2 and from Fig. |SJ we can deduce that for fixed /Edd and r] cm and m < 
10 2 (/co/0.01) 1 / 2 M Q , the duty cycle is larger than unity and, therefore, the background can basically be treated as 
gaussian at all frequencies where a detectable signal is predicted. For to > 10 2 (/co/0.01) 1 / 2 M Q , the frequency where 
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FIG. 7: In the left panel we show the inspiral rate F(z) as a function of redshift for m = 1 — 1O 3 M0. The solid curves represent 
the total plunging rate in the Universe irrespective of frequency, while the dashed curves represent the resolvable fraction with 
maximal SNR > 15 In the right panel we show the total inspiral chirp rate F(/) of Eq. 12911 for the same range of masses, again 
distinguishing between the total rate in the Univers (solid curve) and the resolvable portion (dashed curves). 




FIG. 8: Unresolvable GW spectra for a range of compact object masses m, after subtracting out individual events with 
integrated SNR above a threshold of 15. The solid black curve corresponds to the total time-averaged spectrum, including the 
resolvable portion. 



the duty cycle crosses unity is in a range where T>(f) roughly scales as /~ 8 / 3 (see Fig. EJ), because the inspiral rate 
r(/) varies much more weakly (see Fig. 0), and the coherence time, Eq. I|3U|) . i co h oc J -8 / 3 . As a result, the GW 
background becomes gaussian when 

^/^*10-=(^) 3/8 (^) 3/1 H, (34) 

Note from Fig. [S] that the frequency below which the background becomes unsubtractable, and thus contributes 
to confusion noise, is lower by a factor ~5-10 compared to the frequency Eq. (|34|l below which the noise becomes 
gaussian because the duty factor becomes larger than unity. Figured also shows that the compact object mass above 
which inspirals tend to become resolvable, m > 10 M©, is comparable to the compact object mass above which the 
signal consists of individual events (i.e. the duty cycle is smaller than unity). In other words, there is indeed a range 
of frequencies and compact object masses giving a quasi-continuous signal which may still be resolved into individual 
events given an improved detector. However, when the duty cycle is much greater than unity, even a perfect detector 
will not be able to resolve the individual sources, as is the case for most of the galactic white dwarf population Q- 
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V. 



CONCLUSIONS 



Assuming scaling relations between accretion rates, intrinsic bolometric and X-ray luminosities, and supermassive 
black hole mass on the one hand, and rates and masses of compact object inspirals on the other hand, one can predict 
GW spectra from AGN sources. When part of these scaling relations are constrained observationally, one obtains 
amplitudes detectable by LISA if the fraction of material accreted in the form of compact objects is > 1%. The GW 
signal depends relatively little on the accretion and emission parameters, as long as these are within the observational 
constraints. It increases with the fraction of the observed mass density of supermassive black holes that is due to 
accretion (/Edd), and is thus largest for accretion dominated growth. 

Our scenarios are well within existing observational upper limits. For example, for the highest flux shown in 
Fig. Elone obtains f2 gw (10~ 8 Hz)/i§ < 10~ 12 (/ co /0.01), compared to the msec pulsar timing limit Sl gw (10 — 8 Hz)ft, 2 < 
10~ 8 |3f)|. This class of sources is not important for future missions sensitive around 1 Hz, such as BBO [ifij . 

If LISA establishes a detection or upper limits on fi gw (/), this may be translated into constraints on the parameters 
f co , m, rj gw , fx, /Edd, and r/ cm used in the present parametrization. Since fx, /Edd, and r) em are rather well constrained 
by AGN astronomy, and rj gw is rather well determined by theory, LISA will mainly constrain / co and m. A more 
detailed study where these parameters depend on luminosity and/or different AGN classes may allow for more detailed 
predictions, but is beyond the scope of the present exploratory study. 

For compact object masses m < 10 2 (/co/0.01) 1 / 2 M & , where f co is the fraction of material accreted in form 
of these compact objects, the background is gaussian at nearly all frequencies where a detectable signal is pre- 
dicted. In contrast, for m > 10 2 (/co/0.01) 1 / 2 Mq, the GW background is gaussian only at frequencies / < 
2 x 10~ 3 (/co/0.01) 3 / 8 (10 2 M0/m) 3 / 4 Hz, above which the duty cylce is smaller than unity. The frequency below 
which the GW background becomes confusion noise because the signal to noise ratio becomes too small to subtract 
out individual events is smaller than this gaussian frequency by a factor ~5— 10. At higher frequencies one would ob- 
serve individual coherent signals of typical duration ~ 0.2 (10 2 M©/m) (10 -3 Hz//) 8 / 3 yr which occur at a typical rate 
~ 10 2 (/co/0.01) (10 2 Mq/to) yr -1 . This is typical for supermassive binary inspiral discussed before in the literature. 

While the time-averaged GW background described in this paper will be similar to other EMRI backgrounds 
(e.g. |llj|). the resolvable waveforms should provide information that allows us to distinguish between the circular, 
equatorial inspirals of disk-embedded objects and the highly inclined, eccentric orbits predicted for capture inpirals. 
When combined with knowledge of the electromagnetic properties of active galactic nuclei, LISA should be able to 
constrain primarily the fraction of accretion in form of compact objects and the mass of these compact objects, and 
thus provide fundamental new insights into the problem of star formation in AGN disks. 
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